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As sessile organisms, plants are constantly exposed to a wide spectrunn of stress 
conditions sucli as higin temperature, wliich causes protein misfolding. Misfolded proteins 
are Iniglily toxic and must be efficiently removed to reduce cellular proteotoxic stress 
if restoration of native conformations is unsuccessful. Although selective autophagy is 
known to function in protein quality control by targeting degradation of misfolded and 
potentially toxic proteins, its role and regulation in heat stress responses have not been 
analyzed in crop plants. In the present study, we found that heat stress induced expression 
of autophagy-related (ATG) genes and accumulation of autophagosomes in tomato plants. 
Virus-induced gene silencing (VIGS) of tomato ATG5 and ATG7 genes resulted in increased 
sensitivity of tomato plants to heat stress based on both increased development of 
heat stress symptoms and compromised photosynthetic parameters of heat-stressed 
leaf tissues. Silencing of tomato homologs for the selective autophagy receptor NBR1, 
which targets ubiquitinated protein aggregates, also compromised tomato heat tolerance. 
To better understand the regulation of heat-induced autophagy, we found that silencing 
of tomato ATG5, ATG7, or NBRI compromised heat-induced expression of not only the 
targeted genes but also other autophagy-related genes. Furthermore, we identified two 
tomato genes encoding proteins highly homologous to Arabidopsis WRKY33 transcription 
factor, which has been previously shown to interact physically with an autophagy protein. 
Silencing of tomato WRKY33 genes compromised tomato heat tolerance and reduced 
heat-induced ATG gene expression and autophagosome accumulation. Based on these 
results, we propose that heat-induced autophagy in tomato is subject to cooperative 
regulation by both WRKY33 and ATG proteins and plays a critical role in tomato heat 
tolerance, mostly likely through selective removal of heat-induced protein aggregates. 

Keywords: autophagy, heat tolerance, tomato, WRKY33, NBR1, ATG5, ATG7 



INTRODUCTION 

Autophagy is a highly conserved intracellular degradation system 
in eukaryotes for removal and recycling of cytoplasmic com- 
ponents including damaged proteins and organelles (Klionsky, 
2005). Central to autophagy is the formation of autophagosomes 
resulting from the dynamic membrane reorganization. In yeast, 
more than 30 autophagy-related (ATG) genes have been identified 
and their products often form functional groups that cooper- 
ate to perform the physiologically continuous but mechanisti- 
cally distinct processes of autophagy including the induction of 
autophagy, autophagosome nucleation, elongation, maturation, 
and fusion with vacuoles (He and Klionsky, 2009). Autophagy is 
active at very low levels but is highly inducible in responses to 
stress and extracellular cues (He and Klionsky, 2009). In yeast and 
animal systems, the serine/threonine protein kinase TOR (target 
of rapamycin) functions as a central inhibitor of autophagosome 
formation. In yeast, inhibition of TOR leads to activation of 
ATGl, which can then bind ATG 13 and ATG 17 with increased 
affinities to promote assembly of the ATG1-ATG13-ATG17 scaf- 
fold and initiation of autophagosome formation through recruit- 
ment of multiple ATG proteins (He and Klionsky, 2009). The 



rapid increase in the autophagic flux during the early minutes 
or hours of exposure to stress conditions is mostly mediated 
by post-translational modifications of the core machinery of 
autophagy (He and Klionsky, 2009). A delayed and protracted 
autophagic response, however, also relies on activation of specific 
transcription programs involving stress-responsive transcription 
factors. 

Over the past two decades or so, more than 30 ATG genes 
have been identified in Arabidopsis and other plants including 
tobacco, rice, and maize. Functional analysis of the ATG genes 
has shown that autophagy plays an important role in nutrient 
recycling and utilization in plants (Bassham et al., 2006; Liu 
and Bassham, 2012). Autophagy is also involved in the regu- 
lation of plant senescence, which may be considered a process 
of nutrient redistribution. In addition, autophagy shapes plant 
innate immune responses (Zhou et al, 2014b). Autophagy is 
also induced by a wide spectrum of abiotic stresses including 
oxidative, high salt, osmotic stress and heat conditions (Slavikova 
et al, 2008; Liu et al., 2009). Autophagy-defective mutants or 
transgenic plants are hypersensitive to reactive oxygen species 
(ROS), salt, drought, and heat conditions (Xiong et al., 2007a,b; 
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Liu et al, 2009; Zhou et al., 2013, 2014c). As in other organisms, 
formation of autophagosomes and expression of ATG genes are 
induced by a variety of stresses and environmental cues in plants. 
TOR is also a negative regulator of autophagy in plants (Liu 
and Bassham, 2010). Furthermore, a NADPH oxidase inhibitor 
blocks autophagy induction upon nutrient starvation and salt 
stress, but not during osmotic stress (Liu and Bassham, 2010). 
Thus, ROS may mediate induction of autophagy during some, 
but not aU stress conditions. There is, however, little infor- 
mation available about the transcriptional regulation of plant 
autophagy-associated genes under stress conditions. 

In the present study, we analyze the role and regulation of 
autophagy in heat stress tolerance of tomato plants (Solarium 
lycopersicum) . Heat is an important abiotic stress condition that 
can cause misfolding and denaturation of proteins. As an impor- 
tant protein quality control mechanism, autophagy could play a 
critical role in removal of those misfolded/denatured and poten- 
tially highly toxic proteins or protein aggregates that fail to be 
reestablished for normal protein conformations (Kraft et al, 
2010; lohansen and Lamark, 2011; Shaid et al., 2013). Tomato 
is one of the most important vegetable plants closely related 
to many commercially important plants including potato, egg- 
plant, peppers, tobacco, and petunias and an important model 
plant because it has a number of interesting features such as 
fleshy fruits not shared by other model plants (e.g., Arabidopsis 
and rice) and has many recognized wild species. The tomato 
genome has been sequenced (Consortium, 2012) and its genes 
can be functionally analyzed through a number of complemen- 
tary approaches including virus-induced gene silencing (VIGS) 
(Liu et al, 2002). From the sequenced tomato genome, we have 
identified tomato homologs for ATG5 and ATG7, two important 
autophagy-related proteins that have been subjected to func- 
tional analysis in Arabidopsis (Yoshimoto, 2010; Lai et al, 201 lb; 
Zhou et al, 2013). We have also identified tomato homologs for 
NBRl, a plant selective autophagy receptor (Svenning et al., 201 1; 
Zhou et al, 2013, 2014c), and WRKY33, a transcription factor 
that physically interacts with ATGlSa in Arabidopsis (Lai et al., 
201 lb). Using a variety of molecular approaches including VIGS, 
we analyzed the roles of these genes in heat-induced autophagy 
and heat stress tolerance. These experiments not only support the 
critical role of autophagy in plant heat tolerance but also provide 
new important insights into the regulation of autophagy during 
plant stress responses. 

MATERIALS AND METHODS 

PLANT MATERIALS AND GROWTH CONDITIONS 

Tomato (Solarium lycopersicum L. cv. AQsa Craig) seeds were ger- 
minated in a growth medium fiUed with a mixture of peat and 
vermiculite (7:3, v/v) in trays in a growth chamber. When the 
first true leaf was fully expanded, seedlings were transplanted into 
plastic pots containing the same medium. The growth conditions 
were as follows: light/dark cycle, 22/20° C, and photosynthetic 
photon flux density (PPFD), 600 (imol m"^ s"^. 

QUANTITATIVE RT-PCR (qRT-PCR) 

Total UNA was isolated from tomato leaves using Trizol reagent 
(Sangon Co., Shanghai, China), according to the manufacture's 



recommendations. Genomic DNA was removed with the RNeasy 
Mini Kit (Qiagen Co., Hilden, Germany). 1 |ig RNA was reverse- 
transcribed using the ReverTra Ace qPCR RT Kit (Toyobo Co., 
Osaka, Japan), following the manufacturer's instructions. Gene- 
specific RT-PCR primers were designed based on their cDNA 
sequences (Supplemental Table SI). 

The quantitative real-time PGR was performed using the 
iCycleri QTM real-time PGR detection system (Bio-Rad Co., 
Hercules, CA, USA). Each reaction (25 |xL) consisted 12.5 (jlL of 
SYBR Green PGR Master Mk (Takara Co., Chiga, Japan), 1 |xL 
of diluted cDNA and 0.1|imol forward and reserve primers. 
The PGR cycling conditions and the calculation of relative gene 
expression were as previously described. The tomato ACTIN gene 
was used as internal control as previously described (Zhou et al, 
2014a). 

VIRUS-INDUCED GENE SILENCING (VIGS) 

The tobacco rattle virus (TRV) VIGS constructs for silencing of 
tomato ATG5, ATG7, NBRla, NBRlb, WRKY33a, and WRKY33b 
genes were generated by PGR amplification using gene-specific 
primers (Supplemental Table S2), digested with appropriate 
restriction enzymes and ligated into the same sites of pTRV2. 
The resulting plasmids were transformed into Agrobacterium 
tumefaciens GVdlOl. Agrobacterium-mediatedYims infection was 
performed as previously described (Ekengren et al, 2003). Plants 
were then kept at 22/20° C under 150 [imol m"^ s"' PPFD for 30 
days before they were used for the experiments (Kandoth et al, 
2007). Leaflets in the terminal of the fifth fully expanded leaves, 
which showed 20-30% transcript levels of control plants, were 
used. Each replicate had 12 plants. 

ASSESSMENT OF HEAT TOLERANCE 

To evaluate the role of autophagy in heat stress responses of 
tomato plants, tomato plants at the five-leaf stage were transferred 
to a growth chamber for heat stress treatment (45° G, 400 (imol 
m-^ PPFD, 8 h). 

For determination of electrolyte leakage (EL) caused by high 
temperature, the leaflets in the terminal of the fifth fully expanded 
leaves were measured after heat stress as previous described 
(Huang et al., 2010). Chlorophyll fluorescence was measured 
using an Imaging-PAM Chlorophyll Fluorometer equipped with 
a computer-operated PAM-control unit (IMAG-MAXI; Heinz 
Walz, Effeltrich, Germany). The plants were maintained in the 
dark for more than 30 min before the measurements were per- 
formed. The intensities of the actinic light and saturating light 
were 280 and 2500|xmol mol^^ s^^ PPFD, respectively. The 
maximum quantum yield of PSII (Fv/Fm) was measured and 
calculated as previous described (Zhou et al., 2014a). Three 
replicates for each treatment were used with 12 plants for each 
replicate. 

The light-saturated CO2 assimilation (Asat), stomatal con- 
ductance (Gs), and intracellular CO2 concentration (Ci) were 
determined in the silenced and pTRV plants with an infrared 
gas analyzer-based potable photosynthesis system (LI-6400; Li- 
COR, Lincoln, NE, USA). The air temperature, relative humidity, 
CO2 concentration, and PPFD were maintained at 22° C, 85%, 
380 (imol mol^', and 1000 jxmol m^^ s^', respectively. 
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SEPARATION AND MEASUREMENTS OF TOTAL AND INSOLUBLE 
PROTEINS 

Tomato leaves were collected before and after heat treatment, 
ground in liquid nitrogen and homogenized in a detergent 
containing extraction buffer (100 mMTris/HCl, pH 8.0, 10 mM 
NaCl, 1 mM EDTA, 1% Triton X-100, 0.2% P-mercaptoethanol). 
Soluble and detergent-resistant insoluble proteins were separated 
through low-speed centrifugation and measured as previously 
described (Zhou et al, 2013). 

VISUALIZATION OF INDUCTION OF AUTOPHAGY 

For visualization of autophagosomes, tomato leaves were 
vacuum-infiltrated with 1 \iM of the fluorescence dye 
LysoTracker Green DND-26 (Cell Signaling Technology, 
Danvers, MA, USA) or 500 fiM of mondansylcadaverine (MDC) 
(Sigma- Aldrich, St. Louis, MO, USA). Fluorescence was visu- 
alized using a Zeiss LSM510 UVMeta laser scanning confocal 
microscope (Zeiss Co., Munchen, Germany) and images were 
superimposed using Zeiss LSM510 software. 

RESULTS 

IDENTIFICATION OF TOMATO ATG5. ATG7, AND NBR1 GENES 

To analyze the role of autophagy in tomato heat tolerance, we 
chose first to focus on tomato ATG5 and ATG7 as potential tar- 
gets for gene silencing as their products are required for the 
core process of autophagy and mutants of their Arabidopsis 
homologs, which are single-copy genes, have been widely used 
for functional analysis of autophagy (Yoshimoto, 2010; Lai et al., 
2011b; Zhou et al, 2013). From the sequenced tomato genome, 
we identified two tomato ATG5 genes, ATGSa (Sl02g038380) 
and ATGSb (Sl06g043140). Based on its genomic and full-length 
cDNA sequences, ATGSa has an intron-exon structure similar 
to that of Arabidopsis ATG5 with eight exons and seven introns 
(Supplemental Figure lA). ATGSa encodes a protein containing 
369 amino acid residues and sharing approximately 60% sequence 
identity with Arabidopsis ATG5. The genomic sequence of tomato 
ATGSb is 97% identical to the coding sequence of tomato ATGSa 
(Supplemental Figure IB). The lack of introns in tomato ATGSb 
indicates that the gene was copied from mRNA of tomato ATGSa 
and incorporated into the tomato genome. Further sequence 
analysis revealed multiple mutations in tomato ATGSb, includ- 
ing a G to A transition at nucleotide position 246 that changes a 
tryptophan to a stop codon, resulting in a predicted loss of the 
last 277 amino acids (77%) of the polypeptide if the mutation 
is not removed by alternative splicing (Supplemental Figure 2). 
Thus, tomato ATGSb likely encodes a nonfunctional protein even 
if it is expressed in tomato tissues. Like Arabidopsis, tomato con- 
tains a single ATGZ gene (Slllg068930) that encodes a protein of 
715 amino acid residues. Tomato ATG7 shares approximately 70% 
sequence identity with Arabidopsis ATG7. 

There is a mounting body of evidence that selective autophagy 
of dysfunctional organelles and toxic macromolecules medi- 
ated by selective autophagy receptors play a critical role in 
protein/organelle quality control and in responses to adverse 
environmental and physiological conditions (Kraft et al., 2010; 
Johansen and Lamark, 2011; Floyd et al., 2012; Shaid et al., 
2013). Previously Arabidopsis NBRl, homolog of the mammalian 



autophagy receptors P62 and NBRl, has been analyzed and found 
to play a critical role in plant responses to a spectrum of abiotic 
stresses by targeting stress-induced cytosolic protein aggregates 
(Zhou et al., 2013). Unlike in Arabidopsis, which contains a 
single NBRl gene, tomato contains two NBRl genes, NBRla 
(Sl03gl 12230), and NBRlb (Sl06g071770), with intron-exon 
structures similar to each other and to that of Arabidopsis NBRl. 
Tomato NBRla and NBRlb encodes proteins of 864 and 737 
amino acids, respectively (Supplemental Figure 3), which share 
approximately 50% sequence identify with each other and with 
that of Arabidopsis NBRl. Like Arabidopsis NBRl, both tomato 
NBRla and NBRlb contain two highly conserved ubiquitin- 
associated (UBA) domains and a Wxxl ATG8-interacting motif 
at their respective C-terminus (Supplemental Figure 3). 

HEAT INDUCTION OF AUTOPHAGY-RELATED GENES IN TOMATO 

To determine the role of autophagy in tomato heat tolerance, we 
firstly analyzed the effect of heat stress on expression of tomato 
ATGS, ATG7, and NBRl genes. Tomato seedlings were placed in 
the 22 and 45° C chambers, and the ATG and NBRl gene tran- 
scripts were analyzed by qRT-PCR using total RNA isolated from 
the leaflets of the fifth fully expanded leaves. As shown in Figure 1, 
the levels of tomato ATGS, ATG7 and NBRl transcripts remained 
unchanged at 22°C throughout the 8-h period of the experiments. 
In contrast, transcript levels of the autophagy genes were elevated 
after 2-4 h of heat stress (45°C) and displayed large increases after 
6h of heat stress (Figure 1). After 6-h heat stress, the transcript 
levels of the ATG and NBRl genes started to decline and reduced 
to basal levels by 8 h at 45° C, when the leaves of seedlings started 
to show symptoms of dehydration (Figure 1). Thus, heat stress 
induced the expression of autophagy- related genes. 

SILENCING OF TOMATO ATGS, ATG7, AND NBRl GENES 

To determine directly the role of autophagy in tomato responses 
to heat stress, we used VIGS to assess the impact of down- 
regulated expression of tomato ATGS, ATG7, NBRla, and NBRlb 
on tomato heat tolerance. Gene-specific DNA fragments were 
cloned into the pTRV vector and Agrobacterium cells harbor- 
ing the VIGS vectors were infiltrated into tomato cotyledons. 
We used qRT-PCR to compare the transcript levels for tomato 
ATGS, ATG7, NBRla, and NBRlb in tomato plants infiltrated 
with the pTRV empty vector or infiltrated with the pTRV- ATGS, 
pTRV- ATGZ, plRY-NBRla, or pJRV-NBRlb silencing vector. As 
shown in Figure 2, the basal transcript levels of ATGS, ATG7, 
NBRla, and NBRlb were unchanged in tomato plants after infil- 
tration with the pTRV empty vector but decreased 70-80% in the 
leaves of plants after infiltration with their respective silencing 
vectors (Figure 2). No significant alteration in growth or devel- 
opment was observed upon silencing of the ATGS, ATG7, NBRla, 
or NBRlb gene. 

Abiotic stress including high temperature induces both ATG 
gene expression and formation of autophagosomes. To further 
assess the effect of silencing of tomato ATGS and ATG7 on heat- 
induced autophagy, we used LysoTracker Green dye as a probe 
to detect autolysosome-like structures. The LysoTracker dye has 
been widely used as a probe for detecting autophagic activ- 
ity in a variety of organisms including plants (Otegui et al.. 



www.frontiersin.org 



April 2014 1 Volume 5 | Article 174 | 3 



Zhou et al. 



Autophagy and tomato heat tolerance 




<l> 
o 
c 
CO 

(0 

o 

0) 

c 

CO 



0) <^ 

cc 0 



CD 
O 

I 80 
c 

60 

Q. 



0 2 4 6 
Time (h) 



NBR1a 




0 2 4 6 
Time (li) 



40 



I 20 

TO 

cr 0 



WRKY33a 




2 4 6 
Time 



o 
c 

CC 
"D 

c 

(0 



(0 

DC 0 



ATG7 




-■-22 °C 

0-45 °C 


■ 


0 2 


4 

Time (In) 


6 8 


NBR1b 




-•-22 °C 






-0-45 °C 




■ 


-■ ^ 


^ gr ^ ^ 


0 2 


4 

Time (h) 


6 8 






-«-22 °C 

K>45 °C 




■ 


-■ ■ 








0 2 


4 

Time (h) 


6 8 



FIGURE 1 I Induction of tomato ATG5, ATG7, NBRl, and 
WRKY33 genes by heat stress. Six weeks-old tomato plants 
were placed in the 22 and 45°C growth chambers and total 



RNA was isolated from leaf samples collected at indicated times. 
Transcript levels were determined using real-time qRT-PCR. Error 
bars indicate SE (n = 3). 



2005). Comparative studies with other autophagosome mak- 
ers such as ATG8 have shown that ahhough LysoTracker dyes 
stain acidic organelles, including autophagosomes, up-regulated 
LysoTracker-stained structures are biologically characteristic of 
induced autophagic activity (Phadwal et al, 2012; Chikte et al, 
2014). Under the normal temperature (22°C), we observed 
low numbers of punctate green fluorescent signals in both 
pTRV and gene silenced plants (Figures). After 6-h heat stress, 
however, the numbers of punctate green fluorescent signals 
increased by more than 10 fold in the control plants infil- 
trated with the pTRV empty vector (Figure 3). Importantly, 
in the plants infiltrated with the pTRV-ATGS or pTRV-ATGZ 
silencing vector, there was only a 2-3-fold increase in the num- 
bers of the punctate fluorescence signals after 6-h heat stress 
(Figure 3). 



We also used MDC as a probe for detection of autophagic 
activity in tomato leaves. MDC is an autofluorescent dye that 
stains autophagosomes in mammals and plants (Biederbick et al., 
1995; Munafo and Colombo, 2001; Contento et al, 2005; Liu 
and Bassham, 2010). Under the normal temperature (22°C), 
again, we observed low numbers of punctate fluorescent sig- 
nals in both control (pTRV) and ATG5- or ATGZ-silenced 
plants (Figure 4). After treatment with dithiothreitol (DTT), 
a known inducer of autophagy (Liu et al., 2012), the num- 
bers of punctate fluorescent signals increased by more than 8 
fold in control plants infiltrated with the pTRV empty vector 
(Figure 4). After 6-h heat treatment, the numbers of punc- 
tate fluorescent signals also increased by about 6 fold in con- 
trol plants infiltrated with the pTRV empty vector (Figure 4). 
In the plants infiltrated with the pTRV-ATGS or pTRV-ATGZ 
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silencing vector, the numbers of the punctate fluorescence signals 
after DTT or heat stress were substantially reduced when com- 
pared to those in DTT- or heat-treated control plants (Figure 4). 
These observations confirmed that heat-induced autophagy was 



partially blocked by silencing of the tomato ATG5 and ATG7 
genes. 

COMPROMISED HEAT TOLERANCE OF AUTOPHAGY-SILENCED 
TOMATO PLANTS 

For comparison of heat tolerance of pTRV, pTRV- ATGS, pTRV- 
ATG7, pTRV-NBRla, and pTRV-NBRlb plants, they were placed 
in a 45° C growth chamber for 8h and then moved to room 
temperature for 3-day recovery. For heat-treated pTRV control 
plants, only patches of old leaves displayed symptoms of dehy- 
dration while a majority of the leaves remained green and viable 
after recovery (Figure 5A). On the other hand, a majority of 
fully expanded leaves from the pTRV-ATGS, pJRY-ATG7, pTRV- 
NBRla, and pTRV-NBRlb plants exhibited wilting after the 
recovery (Figure 5A). The more severe symptoms in autophagy- 
suppressed tomato plants after heat stress were confirmed by 
increased EL in the silenced plants relative to that in the pTRV 
control plants (Figure 5B). 

Heat has a harmful effect on various biology processes includ- 
ing photosynthesis. To further investigate responses of ATG5-, 
ATG7-, and NBJ^i -silenced tomato plants to heat stress, we com- 
pared these tomato plants with the pTRV control tomato plants 
for the effects of heat stress on the maximum quantum yield of 
photosystem II (PSII) (Fv/Fm) of leaves immediately after heat 
treatment and light-saturated CO2 assimilation rate (Asat) fol- 
lowing 1-day recovery after heat stress. As shown in Figure 6, 
after 8-h at 45° C, Fv/Fm values for silencing plants were 22- 
37% lower than those of control pTRV plants. Likewise, the Asat 
values for the silenced plants were 35-68% lower than those of 
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FIGURE 2 I TRV-mediated silencing of tomato ATG5, ATG7, NBR1, and 
WRKY33 genes. Transcript levels for each silenced gene in tomato plants 
infiltrated with Agrobactehum cells harboring the empty pTVR vector or the 
corresponding silencing vector were determined by qRT-PCR analysis using 
total RNA isolated from the terminal leaflets of the fifth leaves of 
Agrabscferium-infiltrated tomato plants. Error bars indicate SE (n = 3). 
According to Duncan's multiple range test (P = 0.05), means of lesion areas 
do not differ significantly if they are indicated with the same letter. 



pTRV pTRy- ATG5 pTRV- ATG7 pTRV-WRKYSSa pJR\/-WRKY33b 





FIGURE 3 I Detection of autopfiagic activity using LysoTracker Green. 

Comparison of tomato plants infiltrated with Agrobactehum cells harboring the 
empty pTVR vector or a silencing vector for an indicated gene in terms of 
LysoTracker Green fluorescence signals. The plants were placed into the 22 
and 45°C growth chambers and analyzed after 6-h treatment. Numbers of 



punctate LysoTracker Green fluorescence spots per 10,000 \i,m^ section from 
the cells in the central areas of the terminal leaflets of the fifth leaves were 
indicated below the images. Means and SE were calculated from three 
experiments. According to Duncan's multiple range test (P= 0.05), means of 
lesion areas do not differ significantly if they are indicated with the same letter. 
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FIGURE 4 I MDC-stained autophagosomes in tomato leaves. (A) 

Comparison of tomato plants infiltrated with Agrobacterium cells harboring 
the empty pTVR vector or a silencing vector for an indicated gene in terms 
of MDC fluorescence signals. The plants were placed into the 22 and 45°C 
growth chambers and analyzed after 6-h treatment. DTT (2 mM, 6 h) 
treatment was included as positive control. (B) Numbers of punctate MDC 
fluorescence spots per 10,000 (im^ section were indicated. Means and SE 
were calculated from three experiments. According to Duncan's multiple 
range test (P = 0.05), means of lesion areas do not differ significantly if 
they are indicated with the same letter Bars = 10 (im. 



pTRV control plants when assayed after 1-day recovery following 
heat stress (Figure?). In addition, ATG5-, ATG7-, and NBRl- 
silenced tomato plants had reduced stomatal conductance (Gs) 
and intracellular CO2 concentration (Ci) compared to the unsi- 
lenced pTRV control plants after heat stress (Figure?). Thus, 
photosynthetic efficiency and capacity were more compromised 
by heat stress in the autophagy-suppressed tomato plants than in 
the control plants. 

INCREASED ACCUMULATION OF INSOLUBLE PROTEINS IN 
AUTOPHAGY-SILENCED PLANTS 

Heat stress causes protein misfolding and denaturation, which 
can result in formation of protein aggregates and proteotoxic 
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FIGURE 5 I Functional analysis of tomato ATG5, ATG7, NBR1, and 
WRKY33 in tomato heat tolerance using TRV-mediated gene 
silencing. (A) Tomato plants infiltrated with Agrobacterium cells 
harboring the empty pTVR vector or the silencing pJRV-ATG5, 
pjm-ATG7. pTRV-WBR;, or pjm-WHKY33 vector were placed in a 22 
or 45°C growth chamber for 8h. The pictures of the whole plants 
(upper panel) or the terminal leaflets of the fifth leaves (lower panel) 
were taken after 3-day recovery. (B) Electrolyte leakage (EL) of the 
terminal leaflets of the fifth leaves were determined immediately after 
8h at 22 or 45°C heat treatment. Means and SE were calculated from 
average EL values determined from three experiments with 10 leaves 
per experiment for each genotype. According to Duncan's multiple 
range test (P = 0.05), means of lesion areas do not differ significantly if 
they are indicated with the same letter 
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FIGURE 6 I The effect of heat stress on the efficiency of PSIl 
photochemistry. Fv/Fm images (A) and values (B) of the terminal 
leaflets of the fifth leaves were determined immediately after 8 h at 22 
or 45°C heat treatment. The color code in the images ranged from 0 
(blaci<) to 1.0 (purple). Means and SB were calculated from average 
Fv/Fm values determined from three experiments with 10 leaves per 
experiment for each type of plants. According to Duncan's multiple range 
test (P= 0.05), means of lesion areas do not differ significantly if they 
are indicated with the same letter. 



Stress. To analyze the role of autophagy in protection against 
heat-induced proteotoxic stress, we investigated the accumulation 
of insoluble, detergent-resistant proteins in the pTRV, pTRV- 
ATG5, pTRV-ATG/, pTRY-NBRla, and pTRY-NBRlb plants after 
8 h heat treatment. Total proteins were first isolated and insol- 
uble proteins were separated by low speed centrifugation. As 
shown in Figure 8, the percentages of insoluble to total pro- 
teins were similar in all plants when they were grown at 22° C. 
After 8-h heat stress, insoluble proteins as percentages to total 
proteins increased only by 40% in unsilenced pTRV control 
plants but increased by 138-154% in the pTRV-ATGS, pTRV- 
ATG7, pTRY-NBRla, and pTRY-NBRlb plants (Figures). By 
the end the heat stress, the levels of insoluble proteins in the 
ATG5-, ATG7-, and NBRl -silenced tomato plants were more than 
two times higher than those in the unsilenced control plants 
(Figure 8). 
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FIGURE 7 I The effect of heat stress on the capacity of photosynthesis. 

Light-saturated CO2 assimilation rate (Asat) (A), stomatal conductance (Gs) 
(B) and intracellular CO2 concentration (Ci) (C) were determined following 
1-day recovery after heat stress. Means and SB were calculated from 
average values determined from three experiments with 10 leaves per 
experiment for each type of plants. According to Duncan's multiple range 
test (P = 0.05), means of lesion areas do not differ significantly if they are 
indicated with the same letter. 



IDENTIFICATION AND FUNCTIONAL ANALYSIS OF TOMATO WRKY33 IN 
HEAT TOLERANCE 

Arabidopsis WRKY33 is a transcription factor important for 
plant resistance to necrotrophic fungal pathogens and for plant 
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FIGURE 8 I Accumulation of insoluble protein aggregates under lieat 
stress. Leaf tissues from tomato plants infiltrated with Agrobacterium cells 
harboring the empty pTVR vector or the silencing pJRV-ATGS, pJR\/-ATG7, 
pjm-NBRl, or pJR\/-WRKY33 vector collected at indicated hours (h) under 
45°C for preparation of total, soluble and insoluble proteins as described in 
Materials and Methods. Total proteins in the starting homogenates and 
insoluble proteins in the last pellets were determined the percentages of 
insoluble proteins to total proteins were calculated. 



heat tolerance (Zheng et al, 2006; Li et al., 2011). Arabidopsis 
WRKY33 interacts with Arabidopsis ATGlSa, a critical compo- 
nent of autophagy, and plays a positive role in pathogen-induced 
ATGlSa expression and autophagosome formation (Lai et al, 
2011a). These results suggest that the critical role of WRKY33 in 
plant responses to biotic and abiotic stresses may be at least in 
part mediated through its positive regulation of pathogen/stress- 
induced autophagy. To investigate whether tomato contains 
similar WRKY transcription factor(s) with a critical role in 
heat tolerance and regulation of stress-induced autophagy, we 
searched the sequenced tomato genomes and identified two close 
WRKY33 homologs, WRKY33a (Sl09g0 14990) and WRKY33b 
(Sl06g066370). As shown in Supplemental Figure 4, Arabidopsis 
WRKY33, tomato WRKY33a and WRKY33b all belong to Group 
I WRKY transcription factors containing two WRKY domains 
with highly conserved amino acid sequences. High sequence sim- 
ilarities are also found in the N-terminal domains including 
the highly conserved SP clusters as putative MARK phospho- 
rylation sites and the intervening sequences between the two 
WRKY domains (Supplemental Figure 4). Furthermore, both 
tomato WRKY33a and WRKY33b contain a segment of about 100 
amino acid residues on the C-terminal side of the second WRKY 
domain with substantial sequence homology with ArabidOpsis 
WRKY33, which are absent in other related Group I WRKY tran- 
scription factors such as Arabidopsis WRKY25 and WRKY26 
(Supplemental Figure 4) (Lai et al., 2011a). 

To analyze heat-induced expression of tomato WRKY33 genes, 
we analyzed their transcripts in the tomato seedlings grown 



at 22 or 45° C. As shown in Figure 1, the transcript levels of 
both tomato WRKY33a and WRKY33b remained low through- 
out the 8-h period of the experiments at 22° C. At 45° C, however, 
the transcript levels of tomato WRKY33a and WRKY33b were 
elevated with similar kinetics (Figure 1). Transcripts levels for 
both genes displayed substantial increases after 4-h exposure to 
45°C and peaked after 6-h heat stress (Figure 1). Like those of 
other analyzed autophagy-related genes, the transcript levels for 
both WRKY33a and WRKY33b declined after 6-h heat exposure 
and approached those of control plants by 8-h heat exposure 
(Figure 1). 

To determine directly the roles of tomato WRKY33 genes, 
we used VIGS technology to assess the impact of their down- 
regulated expression on tomato heat tolerance. Tomato WRKY33- 
specific DNA fragments were cloned into the pTRV vector and 
Agrobacterium cells harboring the VIGS vectors were infiltrated 
into tomato leaves. As shown in Figure 2, basal expression of 
WRKY33a or WRKY33b was observed in the tomato plants infil- 
trated with the pTRV empty vector. By contrast, infiltration with 
either p7RY-SlWRKY33a or pTRY-SlWRKY33b silencing vector 
resulted in approximately 5-7-fold reduction in the transcript 
levels for both tomato WRKY33a and WRKY33b (Figure 2). The 
cross silencing likely resulted from the high sequence homology 
between the two genes, which share more than 75% nucleotide 
sequence identify. The tomato plants silenced for WRKY33a and 
WRKY33b were normal in growth and development and dis- 
played no detectable morphological phenotype. 

We analyzed the impact of silencing of WRKY33 genes on 
tomato heat tolerance. Both control and silenced plants were 
placed in a 45° C growth chamber for 8h and then moved to 
room temperature for 3-day recovery. Unlike heat-treated pTRV 
control plants, which had only some patches of old leaves that dis- 
played symptoms of dehydration, a majority of leaves from the 
pTRY-WRKY33a and pTRY-WRKY33b plants exhibited exten- 
sive wilting or even bleaching after the recovery (Figure 5). Thus, 
silencing tomato WRKY33 genes caused increased sensitivity to 
heat stress. Assays of EL, the maximum quantum yield of PSII 
(Fv/Fm), light-saturated CO2 assimilation rate (Asat), stom- 
atal conductance (Gs), and intracellular CO2 concentration (Ci) 
confirmed that silencing of tomato WRKY33a and WRKY33b 
compromised tomato heat tolerance (Figures 5-7). Furthermore, 
silencing of tomato WRKY33a and WRKY33b led to increased 
accumulation of insoluble proteins under heat stress (Figure 8). 

REGULATION OF HEAT-INDUCED AUTOPHAGY BY WRKY33 JKtiD 
AUTOPHAGY PROTEINS 

To investigate whether the critical role of tomato WRKY33 genes 
in heat tolerance is associated with their positive roles in regu- 
lation of heat-induced autophagy, we analyzed whether silenc- 
ing the WRKY33 genes in tomato compromised heat-induced 
autophagosome formation. As shown in Figure 3, while there 
was a more than 10 fold increase in the numbers of punctate 
green fluorescent signals after 6-h heat stress in the pTRV control 
plants, there was only 3-5-fold increase in the pTRV-WRKY33a 
and pTRV-WRKY33b plants. As a result of combined reduction 
of both basal and induced autophagosome formation, the levels of 
autophagosomes in the WRKY33-silenced plants were only about 
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25-30% of those in the control plants after 6-h heat stress. Thus, 
tomato WRKY33 proteins play a positive role in heat-induced 
autophagosome formation. 

We also analyzed the mutual regulation among the silenc- 
ing tomato ATG5, ATG7, NBRl, and WRKY33 genes during 
plant responses to heat stress. For this purpose, we analyzed 
the transcript levels of ATG5, ATG7, NBRl a, NBRlb, WRKY33a, 
and WRKY33b in the silencing plants after 6-h heat stress. 
In pTRV control plants, as expected, the transcript levels for 
ATG5, ATG7, NBRl a, NBRlb, WRKY33a, and WRKY33b were 
elevated under heat stress. However, induction of these genes 
was all reduced not only in the plants harboring their respec- 
tive silencing vectors but also in the plants harboring silencing 
vectors for the other genes (Figure 9). Induction of three heat 



shock proteins {HSP17.6, HSP20, andHSPlOO) in ATG5-, ATG7-, 
NBR1-, or VW^JCY33-silencing plants or pTRV plants was almost 
as strong as that in the pTRV control plants after 6-h heat stress 
(Figure 10). However, induction of a tomato HSP40 gene was sig- 
nificantly compromised by silencing of the autophagy-related or 
WRKY33 genes (Figure 10). These results indicated that ATG and 
WRKY33 proteins have a positive role in heat-induced expression 
of autophagy-related genes. 

DISCUSSION 

In the present study, we analyzed the role of autophagy in 
responses to heat stress in tomato, an important horticultural 
crop. Both the expression ofATGS and ATGZ genes and formation 
of autophagosomes were induced in heat-stressed tomato plants 
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FIGURE 9 I Regulation of heat-induced expression of tomato ATG5, 
ATG7, NBRl, and WRKY33 genes. Tomato plants infiltrated with 
Agrobacterium cells harboring the empty pTVR vector or the silencing 
pTRV-Z^res, pjm-ATG7, pTRV-NBR7, or pjm-WRKY33 vector were 
placed in a 45°C growth chamber and total RNA was isolated from 



leaf samples collected after 6-h heat stress for determination of 
transcript levels of indicated genes by qRT-PCR. Error bars indicate SE 
in = 3). According to Duncan's multiple range test {P = 0.05), means 
of lesion areas do not differ significantly if they are indicated with the 
same letter. 
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FIGURE 10 I Regulation of heat-induced expression of tomato HSP genes. 

Tomato plants infiltrated with Agrobacterium cells harboring the empty pTVR 
vector or the silencing pTm-ATG5, pTRV-AT"G7, pTRV-WBR?, or pTm-WRKY33 
vector were placed in a 45°C growth chamber and total RNA was isolated 



from leaf samples collected after 6-h heat stress for determination of 
transcript levels of indicated genes by qRT-PCR. Error bars indicate SE (n = 3). 
According to Duncan's multiple range test (P = 0.05), means of lesion areas 
do not differ significantly if they are indicated with the same letter. 



(Figures 1, 3, 4). The heat tolerance of autophagy-suppressed 
tomato plants due to silencing of ATG5 and ATG7 genes was 
compromised based on their increased morphological symptoms 
associated with enhanced defects in the efficiency and capacity 
of photosynthesis after heat stress (Figures 5-7). These results 
indicate that autophagy plays an important role in tomato heat 
tolerance. 

Heat stress causes protein misfolding and denaturation. 
Misfolded/denatured proteins are highly toxic due to nonspecific 
binding to a variety of cellular constituent and, therefore, must 
be efficiently removed to prevent proteotoxic stresses (Hightower, 
1991). Extensive studies in yeast and animal organisms have 
revealed that misfolded proteins are recognized by the protein 
quality control system, ubiquitinated by chaperone-dependent 
E3 ubiquitin ligases such as the C-terminus of Hsc70-interacting 
protein (CHIP) and subjected to degradation by the ubiquitin 
proteasome system (UPS) (Kraft et al, 2010; Shaid et al, 2013). 
Very recently we have conducted comprehensive genetic analysis 
of Arabidopsis CHIP E3 ubiquitin ligase and discovered its critical 
role in plant responses to a spectrum of abiotic stresses includ- 
ing heat stress (Zhou et al., 2014c). Previously, it has also been 
reported that Arabidopsis CHIP E3 ubiquitin ligase and Hsc70-4 
mediate plastid-destination precursor degradation through UPS 
when the import of the precursors are blocked in a plastid-import 
mutant (Lee et al, 2009). In the present study, we demonstrated 
that silencing of tomato genes encoding NBRla and NBRlb, 



two close homologs of mammalian ubiquitin-binding autophagy 
receptors P62 and NBRl, also compromised tomato heat toler- 
ance (Figures 5-7). Thus, NBRl -mediated selective autophagy is 
critical in tomato heat tolerance most likely through its activity 
in removing heat-induced misfolded proteins. For degradation 
by UPS, proteins must be unfolded to enter the narrow central 
cavity of its barrel-shaped 20S proteolytic core since the steric 
conditions of a folded protein would not be able to pass through 
the entrance channel. Under heat stress, misfolded, or denatured 
proteins may form protein aggregates that are difficult to dissoci- 
ate or unfold. These protein aggregates are likely to be targeted by 
NBRl -mediated selective autophagy if they fail to be processed 
by UPS. Consistent with this interpretation, compromised heat 
tolerance of NBRl -silenced tomato plants was associated with 
increased accumulation of protein aggregates under heat stress 
(Figures). These results provide further support that UPS and 
NBRl-mediated selective autophagy function in cooperation in 
the removal of misfolded proteins for protection against pro- 
teotoxic stress under adverse environmental and physiological 
conditions (Zhou et al., 2014c). 

Heat stress induced both LysoTracker or MDC-stained 
autolysosome-like structures (Figures 3, 4) and expression of 
ATG and NBRl genes (Figure 1). Interestingly, silencing of an 
ATG or NBRl gene in tomato plants led to down regulation 
of not only the silenced gene but also other sequence-unrelated 
autophagy genes (Figure 9). Since ATG5, ATG7, and NBRl 
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proteins are not known to be regulators of gene transcription, 
their effect on the expression of other genes is probably indi- 
rect and most likely related to induced autophagy under heat 
stress. It is possible that induced autophagy under heat and per- 
haps other stress conditions as well has a positive role in the 
upregulation of autophagy genes. Consistent with the potential 
signaling role of autophagy, silencing of Arabidopsis TOR gene 
leads to not only constitutive formation of autophagosomes but 
also induced expression of some ATG genes (Liu and Bassham, 
2010). We have also previously observed that in autophagy- 
deficient mutants, induction of jasmonate-regulated PDF 1.2 gene 
by Botrytis infection was compromised (Lai et al, 2011b). 
Likewise, in the ATG5- or ATGZ-silenced tomato mutants, induc- 
tion of a gene encoding a HSP40 was significantly reduced 
(Figure 10), indicating that induced autophagy, perhaps through 
formation and turnover of autophagosomes, has a positive role in 
up-regulation of not only autophagy genes but also other genes 
associated with defense and stress responses. We have previously 
shown that in Arabidopsis atgS and atg7 mutants, the protein lev- 
els ofNBRl increased greatly under heat stress (Zhou et al., 2013). 
Other studies have also showed that some of the autophagy- 
related proteins and NBRl are themselves autophagy substrates 
(Svenning et al., 2011) that undergo rapid turnover and must be 
replenished through increased synthesis for sustained autophagy. 
On the other hand, when autophagy is suppressed or blocked as 
in the ATGS- and ATCZ-silenced plants, there would be no need 
for increased transcription of the ATG or NBRl genes since their 
degradation by autophagy is inhibited. 

Arabidopsis WRKY33 transcription factor plays a critical role 
in plant resistance to necrotrophic fungal pathogens and in plant 
tolerance to heat stress (Zheng et al., 2006; Li et al, 2011). 
Autophagy is induced by necrotrophic pathogens, heat, and salt 
stresses and plays an important role in plant responses to these 
biotic and abiotic stresses as well. We have previously shown 
that in the Arabidopsis wrky33 mutants, increased formation of 
autophagosomes was observed in Botrytis-infected lesion areas 
but not in the areas surrounding the lesions found in wild-type 
plants (Lai et al., 2011b). In addition, induction of ATG 18a was 
normal at 1 day post Botrytis infection (dpi) but was severely 
compromised at 2, 3, and 4 dpi in wrky33 (Lai et al., 2011b). 
Thus, WRKY33 is dispensable for early induction of authophagy 
but necessary for sustained induction of autophagy in Botrytis- 
infected plants. Likewise, silencing of tomato WRKY33 com- 
promised heat-induced autophagy gene expression and reduced 
autophagosome formation (Figures 3, 8). Thus, it is likely that the 
critical role of WRKY33 in plant disease resistance and stress tol- 
erance is, at least in part, mediated by its critical role in induction 
of autophagy. 

In Arabidopsis, WRKY33 is subjected to dual-level regulation 
by the mitogen protein kinase 3 and 6 (MPK3/6) cascade (Mao 
et al., 2011). Upon pathogen infection, WRKY33 is phosphory- 
lated by the pathogen/stress-induced MPK3/MPK6 and phospho- 
rylation of WRKY33 is likely to promote the transcription activity 
of WRKY33, which can bind to the W boxes in its promoter and 
turning on its own expression (Mao et al., 2011). Expression of 
Arabidopsis WRKY33 is also induced by abiotic stress conditions 
and by paraquate, which generates ROS in exposed plant cells 



(Zheng et al, 2006). In plants, autophagy is also induced by a vari- 
ety of stresses, including nutrient deprivation, drought, salt stress, 
ROS, and pathogen infection (Liu and Bassham, 2010; Lai et al., 
2011b; Zhou et al., 2013). It has also been shown that NADPH 
oxidase inhibitors block autophagy induction by nutrient starva- 
tion and salt stress, indicating that ROS may also function as a 
signal in induction of autophagy by some environmental stresses 
(Liu and Bassham, 2010). It is tempting to speculate that multi- 
ple stress-initiated pathways may converge to the activation and 
induction of WRKY33 for induction of some autophagy- and 
other stress-related genes. Although autophagy is highly con- 
served in eukaryotic organisms, WRKY transcription factors are 
mostly plant-specific (Zhang and Wang, 2005). A critical role 
of WRKY33 in the regulation of plant autophagy genes would 
strongly indicate that the regulatory mechanisms of autophagy in 
plants have diverged from those in other eukaryotic organisms. 
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Supplemental Figure 1 | Gene structures and coding sequences of tomato 
ATG5a and ATGSb. (A) Gene structures of tomato ATG5a and ATGSb. (B) 
Comparison of tomato ATG5a and ATGSb nucleotide sequences. Identical 
nucleotides are in red. 

Supplemental Figure 2 | Protein sequence comparison of tomato ATGBa 
and ATG5b. Identical amino acid residues between tomato ATG5a and 
ATGSb are in red. 

Supplemental Figure 3 | Protein sequence comparison of tomato NBR1a 
and NBR1b. Identical amino acid residues between tomato BNRIa and 
NBRIb are in red. 

Supplemental Figure 4 | Protein sequence comparison of Arabidopsis 
WRKY33, tomato WRKY33a (SI06g066370), and WRKY33b (SI09g014990). 

Amino acid residues of tomato WRKY33a or WRKY33b identical to those 
of Arabidopsis WRKY33 are in red. The highly conserved WRKYGQK 
sequences and the residues forming the C2H2 zinc fingers are in blue. 
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